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The merging procedure of tree-level matrix elements and the subsequent parton shower as imple- 
mented in the new event generator SHERPA will be validated for the example of VK/Z+jets production 
at the Tevatron. Comparisons with results obtained from other approaches and programs and with 
experimental results clearly show that the merging procedure yields relevant and correct results at 
both the hadron and parton levels. 
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I. INTRODUCTION 

The production of electroweak gauge bosons, e.g. and Z bosons, is one of the most prominent processes at hadron 
colliders. Especially through their leptonic decays they leave a clean signature, namely either one charged lepton 
accompanied by missing energy for W bosons or two oppositely charged leptons for the Z bosons. The combination 
of clear signatures and copious production rates allows a measurement of some of their parameters ^e.g. the W mass 
and width, with a precision comparable with that reached at LEP2 at the Tevatron P,0,0,0,IEIESIEIS EJ > or 
even better at the LHC [TlL IT^ . The same combination, clear signature and large production rate, renders them 
a good candidate process for luminosity measurements, especially at the LHC ^| Q, ITEl llrj| . This holds true in 
particular for TF-bosons, since their production rate is enhanced by roughly an order of magnitude with respect to Z 
production. Furthermore, especially in combination with additional jets, the production of gauge bosons represents a 
serious background to many other interesting processes, leading to multi particle final state topologies. The production 
and decay of pairs of top quarks or of SUSY particles may serve as illustrative examples for such signal processes. The 
special interest in this classic production process is reflected by the fact that it was one of the first to be calculated at 
next-to leading order (NLO) [H III III IS El and next-to-next-to leading order (NNLO) IHH^ in QCD. Recently, 
the first distribution related to these processes, namely the boson rapidity, has been calculated at NNLO |24|. I n 
addition to such fixed-order calculations, programs such as RESBOS [25( have been made available, which resum soft 
gluon effects. Cross sections and distributions for W or Z bosons being produced together with jets can be evaluated 
at the parton level through a number of different computer codes: specialised ones such as VECBOS [2(J, and multi- 
purpose parton-levelgenerators such as COMPHEP [27|, GRACE/GR@PPA [H,^, MADGRAPH/MADEVENT 
|30ll3H. ALPGEN j32|, and AMEGIC [IJ. All of them operate at the tree level, at NLO the program MCFM [HIH 
provides cross sections and distributions for IF/Z+jets for up to 2 jets. 

Apart from such techniques, based on analytical methods, event generators play a major role in the exp erimental 
analysis of collider experiments. In the past years, programs such as PYTHIA |3jj,|37j or HERWIG 38, 39] proved to 
be successful in describing global features of boson production processes, such as the bosons transverse momentum 
or rapidity distribution. Apart from the parton shower, which takes proper care and resums the leading and some 
of the subleading Sudakov logarithms, these programs include the first-order matrix element for the emission of an 
extra parton, implemented through a correction weight on the parton shower. Because of the different ap prox imations 
made for their parton shower, this is realized in different manners inside the two programs, cf. |4jj and |4lll42]|. 
In view of the need for more precise simulations, both in terms of total rates and in the description of exclusive 
final states, two quite orthogonal approaches have been developed recently, which aim at a systematic combination 
of higher-order matrix elements with the parton shower. The first one, called MC@NLO, provides a method to 
consistently match NLO calculations for specific processes with the parton shower. It has been implemented for the 
production of colour-singlet final states, such as W and Z bosons, or pairs of these bosons , or the Higgs boson, and 
for the production of heavy quarks The implementations are available as a code called MC@NLO residing 
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on top of HERWIG. The idea of this approach is to organise the counter-terms necessary to technically cancel real 
and virtual infrared divergencies in such a way that the first emission of the parton shower is recovered. This allows 
the generation of hard kinematics configurations, which can eventually be fed into a parton shower Monte Carlo. The 
alternative approach is to employ matrix elements at the tree level (at leading order) for different jet multiplicities 
and to merge them with the parton shower. The basic idea in this approach is to define a region of jet production, 
i.e. hard emissions, and a region of jet evolution, i.e. soft emissions, divided by a fcj_-type of jet measure [ifil 1471 lisj . 
Leading higher-order effects are added to the various matrix elements by reweighting them with appropriate Sudakov 
form factors. Formal independence at leading logarithmic order of the overall result on the jet measure is achieved by 
vetoing hard emissions inside the parton shower, and suitable starting conditions. This approach was presented for 
the first time [491 for e + e~ collisions; it has been reformulated for dipole cascades in [5(| and extended to hadronic 
collisions in [51j • The method is one of the cornerstones of the new event generator SHERPA > written entirely in 
CH — h, where it has been consistently implemented for nearly arbitrary processes. A case study of this method on the 
basis of PYTHIA and HERWIG, focusing on single-boson prod uction, has been presented in [53. It should be noted 
also that a somewhat related approach has been taken by [53| . There, a systematic mapping and tracing of partons 
stemming from the leading order matrix element through the parton shower takes proper care of leading logarithmic 
effects. 

The goal of this paper is to validate the merging procedure implemented in SHERPA for single-boson production and 
to compare the results with those of other approaches. After a short reminder of the merging procedure and a brief 
introduction to some of the implementation details in Sec.|nJ the focus will shift on results obtained by SHERPA. The 
observables, that will be studied, are inclusive, like the transverse momentum and rapidity distribution of the bosons, 
and more exclusive, like the transverse momentum distribution of additional jets. In a first step, the self-consistency 
of the method will be checked by analysing the dependence of different observables on the separation cut and on the 
maximal number of extra jets provided by the matrix elements, see Sec. lIIII Following this the results of the merging 
method will be contrasted with those of other approaches: on the matrix element level, the jet transverse momentum 
distributions of SHERPAs reweighted matrix elements will be compared with those of a full-fledged NLO calculation 
provided by MCFM; see Sec. lIV Al Then, the results after parton showering and hadronisation will be compared with 
those of other event generators in Sec. lIVBl specifically with those obtained from PYTHIA and MC@NLO. Finally, 
the ability of the method to describe inclusive observables that have been measured, such as the bosons transverse 
momentum, will be exhibited in Sec. II V CI 



II. REALIZATION OF THE MERGING PROCEDURE 



A. Basic concepts 

The key idea of the merging procedure |49l l5l| is to separate the phase space for parton emission into a hard region 
of jet production accounted for by suitable tree-level matrix elements and the softer region of jet evolution covered 
by the parton showers. Then, extra weights are applied on the former and vetoes on the latter, so that the overall 
dependence on the separation cut is minimal. The separation is achieved through a fc^-measure; for hadron collisions 
a longitudinal invariant form is used [tH Esj : two final state particles i and j are defined to be within two different 
jets, if their relative fcj_-measure Qh, given by 

Qij = 2™in{Px,i,P±,j} [cosh(?7 4 - rjj) - cos(<fo - <£,•)] , (1) 

is larger than the jet resolution scale Q 2 ut . In the context of the merging procedure, this scale Q^ ut serves as the 
separation scale between the regimes of jet production and jet evolution. In the above equation, r\ and 4> denote 
the pseudo-rapidities and azimuthal angles of the two particles, respectively. In addition, each jet has to fulfil the 
constraint that 

Ql=Ph (2) 

is larger than the jet resolution scale. Apart from details concerning the merging of two jets into one, variations of 
this jet definition exist. For instance, within the k ± -scheme for Tevatron, Run II, the Qij are eventually rescaled by 
a "cone-like" factor D, and instead of the cosines just the squares of the differences are taken [55j. 
The weight attached to the matrix elements takes into account the terms that would appear in a corresponding 
parton shower evolution. Therefore, a "shower history" is reconstructed by clustering the initial and final state 
particles stemming from the tree-level matrix element according to the fcj_-formalism. This procedure yields nodal 
values, namely the different fc^-measures Q 2 where two jets have been merged into one. These nodal values can be 
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interpreted as the relative transverse momentum describing the jet production or the parton splitting. The four- 
vectors of the mergers are given by the sum of their two offsprings, leading to increasingly massive jets. The first 
ingredients of the ME weight are the strong coupling constants evaluated at the respective nodal values of the various 
parton splittings, divided by the value of the strong coupling constant as used in the evaluation of the matrix element. 
In general, in the matrix element calculation, the jet resolution scale is taken to be the renormalisation scale as well 
as the factorisation scale, guaranteeing that the weight is always smaller than 1. This allows a simple hit-or-miss 
method to yield unweighted events. The other part of the correction weight is provided by Sudakov form factors. At 
next- to leading logarithmic (NLL) accuracy the quark and gluon Sudakov form factors are given by (see |4^) 



A g (Q, Qo) =expl- f dqT q {q, Q) \ , 
I Qo J 



(3) 



A fl (Q,Q )=expJ- J dq[T g (q,Q)+T f (q)]\ , (4) 



Qo 

where Y q , g j are q — » qg, g — > gg and g — > qq branching probabilities 



r, M ) = ^f k «-'|, (5) 



r , ( , iW _4f4(JL"',, ( o) 



7r q \ q 12 

r,(.)-£=f. (7) 

A Sudakov form factor yields the probability for no emission (resolvable at scale Qo) during the evolution from a higher 
scale Q to a lower scale Qo- The ratio of two Sudakov form factors A(Q, Qo)/A(g, Qo) then gives the probability for 
no emission resolvable at a scale Qo during the evolution from Q to q. 

Having reweighted the matrix element, a smooth transition between this and the parton shower region must be 
guaranteed. This is achieved by choosing suitable starting conditions for the shower evolution of the parton ensemble. 
In particular, the starting scale of the shower evolution of a certain parton is not the jet resolution scale Q cu t, but 
rather the production scale of that parton. Then a veto on any emission harder than Q cu t properly separates the 
shower regime from the matrix element region. 



B. The algorithm in general 

The extension of the merging algorithm to hadronic initial states has been proposed in |5lj for the first time. Here, this 
algorithm will be briefly reviewed, with special emphasis on details of its implementation in the SHERPA framework. 
The description of the preferred scale choice for different configurations, details of the treatment of matrix elements 
with the highest jet multiplicity, and the solution to the problem of how to introduce off-shellness to on-shell matrix 
element particles will be considered in the following sections. 
The merging algorithm proceeds as follows: 

1. One process out of all processes under consideration is selected according to the probability 

5>7* 



^ = ~~~~ttY\ ■ (8) 



This choice provides the initial jet rates, subject to an additional Sudakov and coupling weight rejection. For 
instance, a typical selection of processes for W^-boson production would include: 



PP jet jet — > e v e 
pp -» jet jet — > e~v e + jet , 
PP -* jet jet — > e~v e + jet jet , 
pp — * jet jet — > e~v e + jet jet jet . 
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The cross sections <j\ are calculated using the corresponding tree-level matrix elements; the only phase-space 
restriction is given by the fc^-measure. The renormalisation scale [Ir and the factorisation scale fj,p are fixed to 
the cut-off scale Q cu t, with an exception only for the process with the highest number of jets, cf. Sec. Ill Dl 

2. Having chosen a single process, the respective momenta are distributed according to the corresponding differ- 
ential cross section. 

3. The nodal values qi are determined. In doing so a corresponding parton shower history is reconstructed. The 
backward clustering procedure is guided by the fc_L-measure, respecting additional constraints: 

• Unphysical combinations like (qq) and (qq) are ignored. Within the SHERPA framework this is implemented 
by employing the knowledge of the Feynman diagrams contributing to the process under consideration. 
Thus, "unphysical" translates into the non-existence of a corresponding Feynman amplitude. 

• When an outgoing parton of momentum pj is to be clustered with a beam, the fcj_-measure does not provide 
the information as to which beam it has to be clustered. In general the beam with the same sign as its 
longitudinal momentum component is preferred. In addition the new momentum given by p[ — pi — pj 
must exhibit a positive energy in the frame where the initial state shower is performed. 

4. The backward clustering stops with a 2 — > 2 process. The hardest scale of this "core" process has to be found. 
It depends both on the process and its kinematics (cf. Sec. Ill Cl . 

5. The weight is determined, employing the nodal values according to the following rules: 

• For every internal (QCD) line with nodal values qi and qj for its production and its decay, a factor 
A(qi,Q cut )/A(qj,Q cut ) is added. For outgoing lines, a factor A(qi,Q cut ) is added. 

• For every QCD node a factor a s (qi)/a s (Q cu t) is added. 

6. The event is accepted or rejected according to this weight. If the event is rejected, the procedure starts afresh, 
with step^ 

7. The initial and final state showers emerges from the core 2 — > 2 process. Matrix element branchings are included 
as predetermined branchings within the shower. The starting conditions are determined by the clustering 
performed before (in stepEJ), i.e. the evolution of a parton starts at its production scale 1 . The first emission 
from the initial state shower has to take the factorisation scale \xp into account, used during the matrix element 
calculation. 

8. Veto any emission with a k± above the jet resolution scale Qcut- 



C. Special case : W-boson production at hadron colliders 



The algorithm described above and, in particular, the incorporated scale choices, will be illustrated with a few 
examples dealing with VF-boson production: 

The leading order contributions to W~ production are of the Drell-Yan type, i.e. processes of the form 

q <? —> ev e . 

Clustering does not take place, since this is already a 2 — ► 2 process. Furthermore, there is no strong coupling involved, 
and the rejection weight is given by two quark Sudakov factors only: 

W = A q (Q,Q cut )A ql (Q,Q cut ). (9) 

The hard scale Q is fixed by the invariant mass of the fermion pair Q 2 = . Possible configurations resulting from 
the clustering of W + ljet events are exhibited in Fig. ^ The hard 2^2 process either is again a Drell-Yan process 
(Fig. QJ,) or of the type qq' — > gW (Fig-QJ)). The weight in the first case reads: 

W = A q (Q, Q cut ) Ay (Q, Q cut ) A g (Qx, Q cut ) a "^ 1 \ , (10) 



1 Since a virtuality-ordered shower is employed within SHERPA, the virtuality of its predecessor, i.e. its invariant mass, is used. 
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FIG. 1: Two possible cluster configurations of a W+l jet event. The dashed line highlights the hard 2 — » 2 process. 





FIG. 2: Three possible cluster configurations of a W+2 jet event. The dashed line highlights the hard 2 — > 2 process, being 
either of Drell-Yan type (a), a vector boson production (b) or a pure QCD process (c). 



where Q 2 — M e 2 Se and the nodal value Q\ is given by the fc_L-algorithm. For this configuration the gluon jet tends to 
be soft, i.e. Qi preferentially is close to Q C ut- The second configuration differs from the first only by the result of the 
clustering. The transverse momentum of the gluon jet p\ now is of the order of the VF-boson mass or larger. The 
weight looks still the same only the scale definitions are altered. In such a case, the hard scale is now given by 

Q 2 =pl g + M 2 e9e , (11) 

i.e. the transverse mass of the W. Also, the nodal value Q\ has not been determined by the cluster algorithm, 
since it belongs to the (in principle unresolved) core process. A natural choice is the transverse momentum of the 
corresponding jet 

Qi=P±, g - (12) 

These scale definitions guarantee a smooth transition between the two regimes, i.e. from the case where the gluon is 
soft to a case where the gluon is hard. 

More complicated processes involve the production of at least two extra jets. There are many processes contributing 
to this category. Some illustrative examples are displayed in Fig. [21 Cases a) and b) of Fig. [21 are very similar to the 
example with one extra jet only. The corresponding weight reads: 

W = A q (Q,Q cat )A q -,(Q,Q cut )A g (Q u Q cut )A g (Q 2 ,Q cut ) a j^\ ■ ( 13 ) 

ot s {Q cut) a s {Qcut) 



The nodal value Q2 is given by the fcj_-algorithm. The scales Q± and Q are chosen as in the one-jet case. 
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In contrast a new situation arises when a pure QCD process has been chosen as the "core" 2^2 process, see Fig. 
Et) - Since the "core" process is not resolved, there is only one scale available, Q 2 = (2stu)/(s 2 + t 2 + u 2 ) w V\, the 
transverse momentum of the outgoing jets. The correction weight consequently reads: 



W = A g (Q,Q cut ) 



A g (Q,Q cut ) 

Ag(Ql, Qcut) 



A g /(Qi,Q cut ) [A g (Q,Q c 



a s (Q) 
a s (Qcut 



(14) 



The extension to higher multiplicities is straightforward. However, the number of extra jets accounted for by matrix 
elements is limited. This limitation in available MEs enforces a specific treatment of the processes with the highest 
multiplicity. 



D. The highest multiplicity treatment 

In general, the initial cross sections used in step ^ of the merging algorithm above are defined by 

af>= / dxidx 2 dtt fi(xi,/j, F )f 2 {x2,IJ-F)\Mi\ 2 , (15) 



where dil represents the appropriate invariant phase-space element and Mi is the Feynman amplitude for the respec- 
tive process. The choice fip = Q C ut together with the Sudakov factors and the coupling weight leads to a modified 
cross section <n — W af 3 ' . Adding all cross sections with the same number of strong particles yields the cross section 
for production processes accompanied by - exclusively - n jets, 

<&$ = E « • ( 16 ) 

i(njet) 

Of course the number of extra jets that can be considered in this respect is limited by the available matrix elements; 
in SHERPA, this number is usually in the range of three to four. In order to compensate for all the omitted processes 
with more jets, the treatment of processes with the highest number of extra jets differs slightly from the handling of 
lower jet multiplicities. The changes are as follows: 

• the factorisation scale is set dynamically to /ip — <5 m i n , he. to the smallest nodal value as determined by the 
^-algorithm, 

• the resolution scale Q C ut of the Sudakov weights is also replaced by <5 m i n , and 

• the shower veto is applied with Qmin instead of Q C ut ■ 

This guarantees that parton showers attached to matrix elements with the highest number of jets are allowed to 
produce softer jets. In other words: the merging procedure is meant to take into account quantum interference effects 
in jet production at leading order up to a maximal number of jets; any softer jet is left to the parton shower. For the 
configuration shown in Fig. the modified Sudakov and coupling weight reads: 



a s (Qi) a s (Q 2 ) 
a s (Qcut) a s (Q cu t) 



W = A g (Q,Q min )A y (Q,Q min )A fl (Q 1 ,Q min ) "" f ^ L \ ^ \ , (17) 



with lowest scale Q min = Qi- Following this procedure, the sum of all cross sections cfj = VV crf^ for a number of jets 
n can be interpreted as an inclusive cross section 

°t%= E *<■ ( 18 ) 

i(njet) 

i.e. the probability to find at least n jets. Adding all exclusive cross sections for multiplicities lower than a maximal 
multiplicity n max to the inclusive cross section for the highest multiplicity results in a fully inclusive cross section. 
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E. On-shell matrix elements vs. off-shell parton shower kinematics 

One subtle problem when combining matrix elements with parton showers is connected to the question of how to 
translate the matrix element kinematics, determined with on-shell (for light quarks usually massless) momenta, into 
a kinematics suitable for the virtuality-ordered parton shower, i.e. invoking off-shell momenta. Within the SHERPA 
framework, this problem is dealt with by the parton shower in the usual fashion: to begin with, the energy fraction z 
is determined from on-shell kinematics, and afterwards it is reinterpreted. For details on the virtuality-ordered parton 
in SHERPA, the reader is referred to a forthcoming publication |5(tI |. However, if no further emissions are added through 
the parton shower, the partons stemming from the matrix element are left on their mass-shell and the kinematics 
remains unaltered. On the other hand, if the virtuality of one or more partons from the matrix element is increased 
through secondary emissions induced by the shower, the kinematics is modified. Usually, the scales involved in the 
showering are much smaller than the scales prevalent in the matrix elements, which are of the order of or larger than 
Q cu t. Consequently, any manipulation of the kinematics tends to be mild. Nevertheless in some cases changes in the 
kinematics may lead a posteriori to a considerable change of the nodal values from the fc_L-algorithm. For instance, 
the production of W bosons at the Tevatron exhibits a strong asymmetry, which, to a considerable fraction, leads to 
configurations of the initial state, where the emission of jets or extra partons is concentrated on one incoming parton 
only. On rare occasions, such mass effects may alter the smallest nodal value in such a way that it becomes smaller 
than Q cu t. In these cases, the phase-space separation underlying the full merging procedure is violated. Within the 
SHERPA framework, these events are rejected. The corresponding rejection procedure is such that 

• the next event is of the same process as the one being rejected, and 

• the Sudakov weight is only applied to correct the kinematics rather than the rates. 

Therefore, the jet rates determined by the prescription given in the previous section are not altered by the parton 
shower, and the phase-space separation through the /^-algorithm is enforced. 

III. CONSISTENCY CHECKS 

In this section the self-consistency of the results obtained with SHERPA is checked by analysing the dependence of 
different observables on the key parameters of the merging procedure, namely the separation scale Q C ut and the 
highest multiplicity of included matrix elements n max - All plots in this section correspond to W~ boson production 
at the Tevatron, Run II; the parameter settings can be found in the Appendix. If not stated otherwise, the distributions 
shown are inclusive hadron level results, i.e. no cuts have been applied. 

A. Variation of the separation cut Q cu t 

In all figures, the black, solid line represents the total inclusive result as obtained by SHERPA. A vertical dashed line 
indicates the respective separation cut Q cu t, which has been varied between 10 GeV and 50 GeV. To guide the eye, all 
plots also show the same observable as obtained with a separation cut Q cut = 20 GeV, shown as a dashed black curve. 
The coloured lines give the contributions of different multiplicity processes. Note that the separation cut always marks 
the transition between n-jet and n + 1-jet matrix elements. Figs. |3| and 0] show the transverse momentum and the 
rapidity distribution of the W~ boson and the corresponding electron. For the transverse momentum of the W below 
the cut, the distribution is dominated by the LO matrix element with no extra jet, i.e. the transverse momentum is 
generated by the initial state parton shower only. Around the cut, a small dip is visible in Fig. EI T ne p± distribution 
of the electron, in contrast, is hardly altered. The rapidity distributions in Fig. 01 exhibit the asymmetry, which has 
been anticipated when considering merely the negatively charged W boson. The shape of these distributions is very 
stable under a variation of the separation cut. In all observables a small increase of the total cross section of a few 
percent when changing <3 cut from 10 GeV to 50 GeV is visible. This underlines the fact that the dependence on the 
separation cut is weak. 

Differential jet rates with respect to the fc^-algorithm are interesting observables, since they basically exhibit the 
distributions of nodal values using the cluster algorithm. For simplicity the Run II fc^-algorithm has been used with 
D = 1 for the analysis. Differential jet rates are of special interest, since the nodal values are very close to the measure 
used to separate matrix elements from parton shower emissions. Some minor problems with respect to the separation 
should immediately manifest themselves in these distributions. In Fig. [S] the 1 — > 0, 2 — > 1 and 3^2 differential jet 
rates are shown. Within the given approximation the independence is satisfactory. 
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FIG. 3: p±(W~) and p ± (e~) for Q cut = 10 GeV, 30 GeV and 50 GeV in comparison with Q cut = 20 GeV. 

B. Variation of the maximal jet multiplicity n max 

For very inclusive observables such as transverse momentum and rapidity of the W boson, it is usually sufficient to 
include the matrix element with only one extra jet in order to obtain a reliable prediction. Consequently, the inclusion 
of matrix elements with more than one extra jet in the simulation should not significantly change the result. This can 
be used as another consistency check. Figs. Eland impressively picture the dependence on the maximal jet number 
in the matrix elements included. They show that the treatment of the highest multiplicity (cf. Sec. Ill D^i completely 
compensates for the missing matrix elements, whereas the contribution of the lowest multiplicity is not altered. 

C. Matrix element, parton shower and hadronisation 

In addition to the self-consistency of the algorithm tested so far at the hadron level, it is worth while to check that 
the parton shower and hadronisation do not induce significant changes with respect to the initial reweighted matrix 
element in high-p^ regions. Fig. [S] proves that the predictions of SHERPA, e.g. the p± distribution of the hardest jet 
in W production, are remarkably stable in the region of matrix element dominance. 

D. Variation of factorisation and renormalisation scale 

Finally, the sensitivity of the previous results with respect to changes in the renormalisation and factorisation scale 
are examined. In the following all scales occurring in the event generation, both on the matrix element and at the 
parton shower level, are multiplied by constant factors, ranging from 0.5 up to 5. It is clear that the total cross section 
changes with changing scales: starting with 930 pb for the default scale choice, 887 pb (959 pb) are obtained when 
using a scale factor of 0.5 (2). The shape of the p± distributions of the jets, however, experiences only mild changes. 
This is greatly exemplified by the left panel of Fig. |5J where the p± spectrum of the hardest jet is displayed. In the 
right panel of Fig. the result of the left panel is broken down for two different scale prefactors, 1 and 5, to the 
different contributions. Clearly, at the individual level different jet multiplicities differ also in their shapes; in their 
interplay, however, these effects cancel in terms of the overall shape. 
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FIG. 4: n(W-) and r){e~) for Q cut = 10 GeV, 30 GeV and 50 GeV in comparison with Q cut = 20 GeV. 



IV. SHERPA VS. DATA AND OTHER MCS 

In order to study the impact of the merging prescription, the predictions obtained with SHERPA may also be compared 
with other approaches. In a first step, in Sec. IIV Al the transverse momentum distribution of the first and second 
hardest jets in exclusive and inclusive boson-plus-jet production from SHERPA are confronted with the NLO QCD 
predictions of the parton level generator MCFM [53, |3^| . In Sec. IIV Bl the full event generators MC@NLO [t^ and 
PYTHIA 37] arc used to investigate the capabilities of SHERPA when studying W/Z+jet production at the hadron level. 
Finally Sec. IIV CI contains a comparison of the predictions made for the bosons transverse momentum distribution 
with those measured by the DO and CDF collaborations at the Tevatron, Run I. 



SHERPA vs. MCFM 



In order to compare the SHERPA predictions for W/Z+ljet and W/Z +2jet production, a two-step procedure is chosen. 
In a first step the Sudakov and a s reweighted matrix elements are compared with exclusive NLO results obtained 
with MCFM. In the case of the next-to-leading order calculation, the exclusiveness of the final states boils down to 
a constraint on the phase space for the real parton emission. The exclusive SHERPA results consist of appropriate 
leading order matrix elements with scales set according to the fcj_-clustering algorithm and made exclusive by suitable 
Sudakov form factors, cf. Sec. Ill Bl In a second step, the jet spectra for inclusive production processes are compared. 
For the next-to-leading order calculation, this time the phase space for real parton emission is not restricted and the 
SHERPA predictions are obtained from a fully inclusive sample, using matrix elements with up to two extra jets and 
the parton showers attached. If not stated otherwise, all results have been obtained using the input parameters and 
phase-space cuts summarised in the Appendix. Jets are found using the Run II k± -clustering algorithm defined in 
55] with a pseudo-cone size of D = 0.7 and a minimal p± of 15 GeV. 



10 




log(Q| / GeV) \og{Q, / GeV) \og(Q, / GeV) 

FIG. 5: Differential jet rates for the 1 — > 0, 2 — > 1 and 3^2 transition (top to bottom), for 

Qcut = 10 GeV, 30 GeV, and 50 GeV (from left to right). In each plot, the results are compared with those for Q C ut = 20 GeV. 

1 . Exclusive jet p± spectra 

In Fig.^Jthe jet p± distribution for the exclusive production of W+ljet and Z+ljet are shown. In both figures, the 
SHERPA prediction is compared with the exclusive NLO result obtained with MCFM and with the naive LO prediction, 
which is the same for the two programs. For the fixed-order NLO and LO result, the renormalisation and factorisation 
scales have been set to [ir — /if — 80.419 GeV = Mw- All distributions have been normalised to the corresponding 
total cross section. This allows for a direct comparison of the distributions shape. As stated above, the SHERPA results 
stem from Sudakov and a s reweighted W+ljet or Z+ljet LO matrix elements. The change between the naive leading 
order and the next-to-leading order distribution is significant. At next-to-leading order the distributions become much 
softer. For a high-pj^ jet it is much more likely to emit a parton that fulfils the jet criteria and therefore removes 
the event from the exclusive sample. The SHERPA predictions show the same feature. The inclusion of Sudakov form 
factors and the scale setting according to the merging prescription improves the LO prediction, resulting in a rather 
good agreement with the next-to-leading order result. 

In the high-p^ tail, however, the NLO calculations from MCFM tend to be a bit below the SHERPA results. The 
reason is simply connected to the fact that relevant scales in the high-pj^ tail are much larger than the default choice 
of hr = fiF = 80.419 GeV. In order to highlight this, Fig. ^] contains the jet p± distribution in Z+ljet events. In 
this plot, the renormalisation and factorisation scales have been chosen to be fiR = fip = 160.838 GeV = 2Mw- 
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FIG. 6: p±(W ) for Q cu t = 15 GeV and different maximal numbers of ME jets included. 




FIG. 7: r/(W ) for Q cu t = 15 GeV and different maximal numbers of ME jets included. 



Changing the scale in this manner indeed has quite a small impact on the total cross section at NLO, but the tail of 
the distribution becomes considerably enhanced. With the above choice of /i/? and \ip the agreement of NLO and the 
SHERPA result is impressive. 

The p± distribution of the first and second jets in I / F+2jet and Z+2jet production are presented in Fig. Again, 
the next-to-leading oder distributions are softer than the leading order ones, for the same reason as for the ljet case. 
In addition, at low-p^ the leading order result is smaller than the ncxt-to-lcading order one. Taken together, the 
curves have a significantly different shape over the whole interval. This situation clearly forbids the use of constant 
if-factors in order to match the leading order with the next-to-leading order result. Nevertheless, as before, the 
SHERPA prediction reproduces to a very good approximation the shape of the NLO result delivered by MCFM. Fig. 
1131 shows that, similar to the Z+ljet case for VF+2jet in the high-p^ tail, the situation is even better using higher 
renormalisation and factorisation scales (e.g. = /.if = 160.838 GeV) in the NLO calculation. 



2. Inclusive jet p± spectra 

NLO results for inclusive boson plus jet(s) production obtained with MCFM are compared with fully inclusive samples 
generated with SHERPA. There, the matrix elements for W/Z+0,l,2jet production have been used including the highest 
multiplicity treatment for the W/Z +2jet case. The Sudakov and a s reweighted matrix elements have now been 
combined with the initial and final state parton showers. The hadronisation phase for the SHERPA events has been 
discarded. As for the exclusive case the naive leading order prediction is given by the corresponding leading order 
matrix clement that is identical to the one in Figs. Un|and[T21 For the NLO prediction again the renormalisation and 
factorisation scales have been chosen to coincide, namely [1r = fj,p = 80.419 GeV. 

In Fig. ^] the p± spectra for the hardest jet in inclusive W/Z+ljet production are shown. Compared with the 
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FIG. 8: The p± of the hardest jet in inclusive W production at the Tevatron, Run II. The solid line indicates the distribution 
as delivered by the matrix elements; the dashed line is obtained after parton shower evolution; and the dotted line gives the 
final result after hadronisation. Here Q cu t = 15 GeV. 




FIG. 9: The p± of the hardest jet in inclusive W production at the Tevatron, Run II, and its dependence on different choices 
in the renormalisation and factorisation scale. In the left plot, the solid line indicates the default hadron level result; the green 
dashed and the red dotted line are obtained by multiplying all scales in the coupling constants and the PDFs in both the matrix 
elements and the parton shower by a factor of 0.5 or 2, respectively. In the right plot, the same observable is investigated 
in more detail with the default scales (solid lines) and with all scales multiplied by a factor of 5. In addition the different 
contributions of different jet multiplicities are shown. 



exclusive predictions, the high-p_i_ tail is filled again and, hence, the differences between the NLO calculations and the 
LO ones appear to be smaller. For both cases the SHERPA result and the NLO calculation are in good agreement. 
In Fig. ^3 the p± spectra for the first and second hardest jets in inclusive W/Z+2jet production are presented. 
Considering the scale dependence of the next-to-leading order result in the high-p^ region, as already studied in Fig. 
1131 for the exclusive result, the curves are in pretty good agreement. 

Altogether, the merging procedure in SHERPA, including the scale-setting prescription of the approach and the Sudakov 
reweighting of the LO matrix elements, proves to lead to a significantly improved leading order prediction. Seemingly, 
it takes proper care of the most relevant contributions of higher order corrections. Although it should be stressed that 
the rate predicted by SHERPA is still a leading order value only, a constant .ff-factor is sufficient to recover excellent 
agreement with a full next-to-leading order calculation for the distributions considered. Furthermore, by looking at 
the inclusive spectra it is obvious that this statement still holds true after the inclusion of parton showers and the 
merging of exclusive matrix elements of different jet multiplicities. 
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FIG. 10: Jet p± distribution of exclusive W + ljet (left) or Z + ljet (right) events at the Tevatron, Run II. 
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FIG. 11: Jet p± distribution of Z + ljet events at the Tevatron where for the NLO and LO calculation the renormalisation 
and factorisation scales have been chosen to be (Xr = plf = 160.838 GeV. 
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FIG. 12: The p± distribution of the first and second jets in exclusive W + 2jet (left) and in exclusive Z + 2jet (right) events 
at the Tevatron, Run II. 
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FIG. 13: The p± distribution of the first and second jets in exclusive W + 2jet events at the Tevatron where for the NLO and 
LO calculation the renormalisation and factorisation scales have been chosen to /j,r — pf — 160.838 GeV. 
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FIG. 14: The p± distribution of the hardest jet for inclusive W + ljet (left) and for inclusive Z + ljet (right) production at 
the Tevatron, Run II. 
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FIG. 15: The p± distribution of the hardest two jets for inclusive W + 2jet (left) and for inclusive Z + 2jet (right) production 
at the Tevatron, Run II. 



15 




FIG. 16: The p± (left) and r\ (right) distribution of the W boson for inclusive production at the Tevatron, Run II. Plotted 
are the results as obtained with the generators PYTHIA (green), MC@NLO (red) and SHERPA (black). 

B. SHERPA vs. MC@NLO and PYTHIA 

In a next step, results obtained with SHERPA arc compared with those obtained from two other event generators, namely 
PYTHIA and MC@NLO. In both, standard settings have been used for inclusive W production at the Tevatron and 
the underlying event has been switched off. The corresponding process number in PYTHIA is MSEL=12, the relevant 
MC@NLO process is IPR0C=-1471. Inclusive quantities like, for instance, the p± and r\ distributions of the W are in 
good agreement, see Fig. 1161 and only in the high-p^ tail of the distribution some small deviations become visible. 
However, more exclusive quantities such as the pj_-distributions of the first three jets show differences that increase 
with the increasing order of the jet. This can clearly be seen from Fig. El The predictions for the hardest jet start 
to disagree with a factor of roughly 2 at )et-p±s of the order of 100 GeV, reaching up to nearly an order of magnitude 
at p± around 200 GeV. This trend is greatly enhanced for the second and third jets, where discrepancies are of the 
order of one magnitude for the second jet at p± « 100 GeV or even higher for the third jet. 

These discrepancies, however, were to be expected since the other two programs do not include any higher order 
correction beyond first order in the strong coupling constant. In the case of MC@NLO, predictions have been 
compared with those obtained from MCFM; after a careful calibration of input parameters such as CKM elements 
etc., inside the code both programs coincided in all observables tested |57| . Therefore, differences in the p± distribution 
of the hardest jet have to be attributed to a combination of distinct parameter settings and of differences in parton 
showering and hadronisation. The latter type of difference should be taken as some kind of theoretical uncertainty. 
For higher jet configurations, however, the remaining discrepancies are due to different physics input. 




FIG. 17: Jet p± distribution of the three hardest jets in inclusive W production at the Tevatron, Run II. Compared are the 
hadron level results of SHERPA (black), PYTHIA (green) and MC@NLO (red) after 2.5 million events. 
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C. SHERPA vs. data 



Having compared the SHERPA predictions for the case of the jet transverse momentum distributions in exclusive and 
inclusive W/ 'Z+ljet and W/Z+2jet production against other Monte Carlo programs, a comparison with experimental 
data provides an ultimate test of SHERPA's ability to describe such processes. Unfortunately, so far only the inclusive 
W- and Z-boson transverse momentum distribution measured in Tevatron, Run I, have been published, which allows 
for an overall check only. In both cases, matrix elements with up to four (W) or three (Z) extra jets have been taken 
into account - as indicated by the different colours - to generate the SHERPA sample. The black line represents the 
sum of all contributions. For this sample the required separation cut has been chosen to Q cut = 20 GeV. 
In Fig. ^3 the (inclusive) p± distribution of the W is compared with data from DO, taken at Run I of the Tevatron 
[Hsj . The agreement with data is excellent. It can be recognised that approaching the merging scale from below, the 
W+Ojet contribution steeply falls and the distribution for larger momenta is mainly covered by the VF+ljet part, as 
expected. In order to match the measured distribution, the SHERPA result has been multiplied by a constant A-factor 
of 1.25. 

Similarly, in Fig. 1191 the (inclusive) p± distribution of the Z is compared with data, this time taken by CDF at Run 
I of Tevatron 59] . Again the overall agreement is excellent. This time the result has been multiplied by a constant 
A"- factor of 1.6 to match the data. The result is perfectly smooth around the merging scale of Q C ut = 20 GeV. This 
is especially highlighted in the left plot of Fig. ^jj which concentrates on the low momentum region. It is interesting 
to note that the description of the data for momenta smaller than the merging scale is almost only covered by the 
Z+Ojet contribution and is therefore very sensitive to the details of the parton showers and the treatment of beam 
remnants. A parameter of specific impact on the very low momentum region therefore is the primordial (or intrinsic) 
k± used for the interacting partons. This is modelled through a Gaussian distribution with a central value of 0.8 
GeV. Nevertheless, the shower performance of SHERPA has not been especially tuned; the low momentum behaviour 
may therefore still be improved once a detailed parameter tune is available. 




60 80 100 120 140 160 180 200 



P iW / GeV 



FIG. 18: The p± distribution of the VF-boson in comparison with data from DO at the Tevatron, Run I |5fil |. The total result 
is indicated by the black line. The coloured lines show the contributions of the different multiplicity processes. Here matrix 
elements with up to four extra jets have been considered. The applied separation cut is Q C ut = 20 GeV. 
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FIG. 19: The p± distribution of the Z-boson in comparison with data from CDF at the Tevatron, Run I [5flj . The total result 
is indicated by the black line. The coloured lines show the contributions of the different multiplicity processes. The applied 
separation cut is Q C ut = 20 GeV. The right plot focuses on the low momentum region of the left one. 

V. CONCLUSION 

In this work, predictions for single-boson production processes at the Tevatron have been obtained from the new event 
generator SHERPA and compared with the results from other programs and with data. In all cases an encouraging 
agreement of results has been found. Especially for the description of exclusive multi particle final states, SHERPA 
proved its unique value as the simulation tool of choice, limited only by the maximal number of external particles 
covered by its intrinsic matrix elements. In practical applications the choice of this input should be guided by the 
question in consideration. For instance, for the simulation of W + 4jet background to top pair production, a highest 
jet multiplicity of 4 is advisable in order to include most of the quantum interferences. 

Having validated the versatility of SHERPA on one of the most important processes at hadron colliders, its abilities will 
be further tested in the near future by considering other important processes, such as the production of boson pairs, 
the Higgs boson, heavy quarks, or multijet final states. 
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APPENDIX A: INPUT PARAMETERS AND PHASE-SPACE CUTS 



The PDF set used for all analyses is cteq61 The value of a s is chosen according to the value taken for the PDF, 
namely 0.118. For the running of the strong coupling the corresponding two-loop equation is used. Jets or initial 
partons are restricted to the light flavour sector, namely g,u, d, s, c. In fact these flavours are taken to be massless 
and the Yukawa couplings of the quarks are neglected throughout the entire analysis. 



a. SM input parameters 



The SM parameters are given in the scheme: 



m w = 80.419 GeV , T w = 2.06 GeV, 
m z =91.188 GeV, T z = 2.49 GeV, 
G M = 1.16639 x 10~ 5 GeV~ 2 , 
sin 2 W = 1- m^/m| = 0.2222, 

a s = 0.118. (Al) 
The electromagnetic coupling is derived from the Fermi constant G M according to 

The constant widths of the electroweak gauge bosons are introduced via the fixed-width scheme. CKM mixing of the 
quark generations is neglected. 



b. Cuts and jet criteria 

For all jet analyses the Run II fc^-clustering algorithm defined in [55| is used. The parameter of this jet algorithm is 
a pseudo-cone of size D given below for the Tevatron analysis. For the charged leptons the following cuts are applied: 

plepton > 2Q QeV ^ i^epto^ < ^ m ll > lg Ge y ( A3 ) 

For the case of W production an additional cut on missing transverse momentum according to the neutrino has been 
required, namely 

For the jet definition a pseudo-cone size of D = 0.7 has been used in addition to cuts on pseudo-rapidity and transverse 
momentum: 

> 15 GeV, |ry jct | < 2. (A5) 
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